Unconventional decay law for excited states in closed many-body systems 
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We study the time evolution of an initially excited many-body state in a finite system of interacting 
Fermi-particles in the situation when the interaction gives rise to the "chaotic" structure of com- 
pound states. This situation is generic for highly excited many-particle states in quantum systems, 
such as heavy nuclei, complex atoms, quantum dots, spin systems, and quantum computers. For a 
strong interaction the leading term for the return probability W^(i) has the form W(t) — exp(— Aljt'^) 
with A|; as the variance of the strength function. The conventional exponential linear dependence 
Wit) — Cexp(— Ff) formally arises for a very large time. However, the prefactor C turns out to be 
exponentially large, thus resulting in a strong difference from the conventional estimate for W{t). 
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It is known that highly excited states can be treated 
as "chaotic" ones in many-body systems, such as com- 
plex atoms JH , multicharged ions |^ , nuclei [|| and spin 
systems |^,^ This happens due to a very high density 
of many-particle states, which strongly increases with an 
increase of energy. For example, in the case of n Fermi- 
particles occupying the finite number m of "orbitals" 
(single-particle states) , the total number N of many-body 
states grows exponentially fast with an increase of the 
number of particles, N = ml/n\(m—n)\ ~ exp(con). Cor- 
respondingly, the density p/ of those many-body states 
which are directly coupled by a two-body interaction, also 
grows very fast Therefore, even a relatively weak 
interaction between the particles can lead to a strong 
mixing between unperturbed many-body states ("basis 
states"). As a result, an exact (perturbed) eigenstate is 
represented by a chaotic superposition of a large number 
of components of basis states. 

The number of principal basis components in such 
chaotic eigenstates can be estimated as Npc ~ T/D where 
r is the spreading width of a typical component, that can 
be estimated using the Fermi golden rule, and D~-^{E) 
is the total density of many-body states. In the case of a 
quantum computer, the interval between multiqubit en- 
ergy levels D (X 1/iV is extremely small, and practically 
it is impossible to resolve these levels. Moreover, both 
the temperature and finite time of computer operations 
lead to an energy uncertainty 6E 3> £>. A similar situa- 
tion occurs for an electron which enters a many-electron 
quantum dot. In these cases the analysis of stationary 
chaotic eigenstates is not an adequate to real physical 
problems and one needs to consider the time evolution 
of wave functions. In this Letter we extend the quantum 



chaos approach to the problem of time evolution of an 
initially excited basis state. 

Exact many-body eigenstates |fc) of the Hamiltonian 
H = Hq + V oi interacting Fermi-particles can be ex- 
pressed in terms of simple shell-model basis states \f) of 

|/)=4...a+ |0). (1) 

/ 

Here |0) is the ground state, af is the creation operator 
and C'j'^^ are components of an exact eigenstate in the 
unperturbed basis. 

In application to quantum computer models the Hamil- 
tonian Hq describes a number of non-interacting quhits 
(two- level systems), and V stands for the interqubit in- 
teraction needed for a quantum computation ( we assume 
time- independent V). In this case the basis state |/ > 
is a product of single qubit states, is the spin-raising 
operator (if the ground state |0) corresponds to spins 
"down"), and chaotic eigenstates \k) are formed by the 
residual interaction V . 

Below we consider the time evolution of the system, 
assuming that initially {t — 0) the system is in a specific 
basis state \i) (in the state with certain spins "up" for 
a quantum computer). This state can be expressed as a 
sum over exact eigenstates, 

k 

therefore, the time-dependent wave function reads as 

= E ^i'^^? I/) exp(-^£;(^■)^). (3) 
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Here E'^''^ are the eigenvalues corresponding to the eigen- 
states \k >. The sum is taken over the eigenstates \k > 
and basis states |/ > (in what follows, we put h — 1). 

The probability Wi = \A,\^ = | p to find the 

system in the state |i > is determined by the amplitude 



J P^{E)cxp{-iEt)dE. (4) 



Here we replaced the summation over a large number 
of the eigenstates by the integration over their energies 
E = E^''\ and introduced the strength function (SF) 
Pi{E) which is also known in the literature as the local 
spectral density of states, 



(5) 



Here p{E) is the density of states of the total Hamilto- 
nian H, and the average is performed over a number of 
states with close energies. 

In chaotic systems the strength function Pi{E) is 
known to have the Breit-Wigner form for a relatively 
weak interaction, and is close to the Gaussian for a strong 
interaction [|l],D. Recently, the following approximate 
general expression has been analytically found 
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space. It is easy to find the evolution of Wi(t) on a small 
time scale. Let us subtract the energy Ei = Ha of the 
initial state in the exponent and make a second order 
expansion in E — Ei in Eq.lQ). This gives the following 
result. 



and 



A, = exp(-^^,^) (1 - Alt^/2) 



(8) 



(9) 



Here the width A^; of the SF is determined through 

r2 

model is Hctl, 



= E which for the TBRI 
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VqU^ti — l){m — n){m — n + 3), 



(10) 



with Vq standing for the variance of the off-diagonal ele- 
ments (A|; for a quantum computer model can be found 
in Ref. g). 

Note that for a strong residual interaction, Fq > A^, 
the time dependence (^ is also correct for a longer time 
Indeed, both the strength function and density of 
states in this limit are described by the Gaussian func- 
tions with the variance = A|, (see details in ||] - Q), 
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r.,{E)o,27T\V,f\^Pf{E). (7) 

which is derived by making use of the approach described 
in Ref. [^j. Here ^i{E) is some function of the total 
energy. Si is the correction to the unperturbed energy 
level Ei due to the residual interaction Vif, and Pf{E) 
is the density of the basis states |/) directly connected 
with a given state |i) by the matrix element Vif. The 
above result has been derived for the so-called Two-Body 
Random Interaction (TBRI) model ||^ which describes n 
interacting Fermi-particles distributed over m orbitals, 
with an assumption that two-body matrix elements are 
completely random. 

It is shown that for a large number of particles the 
function ViiE) has the Gaussian form with the variance 
which depends on the model parameters. In the case of a 
relatively small (but non-perturbative) interaction (when 
Fj = 2tt pfVi < 2AB with = {V^f)), the function 
T{E) is very broad (i.e. it does not change significantly 
within the energy intervals F and A^;) and can be 
treated as constant, F(£') ~ Fq. In the other limit case 
of a strong interaction, Fo > As, the dependence F(i?) 
in ^ is the leading one. 

The knowledge of the strength function allows one to 
describe the dynamics of wave packets in the energy 



N f E^\ 

Therefore, at the center of the energy spectrum, Ec = 0, 
Eq. (^ results in the Gaussian time dependence for Ai (t) 
and W^{t), 

A, = exp {-A%t^/2) , (13) 

VF,(t)~exp(-A|t2). (14) 

Now let us consider large times. In this limit the result 
can be obtained by evaluation of the integral in Eq. (^) in 
the complex iJ-plane. Specifically, one should close the 
contour of integration in the bottom part of the complex 
plane (Im E < 0), in order to provide a vanishing con- 
tribution at infinity. Then, the large time limit is given 
by the pole of the strength function (|6|), closest to the 
real i?-axis. If F^ and St in Eq.(^) do not depend on E, 
the integration gives the conventional exponential decay 
Wi = exp(— Fi) 1^. However, the energy dependence of 
F is necessary to provide the finite second moment A|; 
of the strength function. If F < A^, the closest pole is 
given by F = —2ImEp, where Ep is the solution of the 
equation Ep = Ei + Si{Ep) — iT{Ep)/2 with a minimal 



imaginary part. If F ^ A^, then we have F w F. As a 
result, we obtain an exponential dependence for large t, 



W^{t) = Cexp(-n), 



(15) 



with some constant C. 

In Ref . jl^ simple extrapolation formula for Wi {t) has 
been suggested, that interpolates (for F < A^) between 
the small (0) and large (p^ time dependencies. 



W^it) = exp 
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(16) 



The transition between the Gaussian regime and simple 
exponential decay occurs near the time ~ T/ where 

2Al 



A 2 y-2 



Ttc- This gives an estimate, C exp 



(2^)' 



for 



the constant C. Indeed, for t < tc the quadratic ex- 
ponential decay, exp(— Al^t^), is slower than the linear 
one, exp(— Fi). The matching of these two dependen- 
cies would naturally require the above expression for C. 
Thus, the constant C can be large if F > Ae- The tran- 
sition from one regime of the time dependence of Wi{t) 
to another is schematically shown in Fig.l. 



It can be used, in conjunction with Eq.(^), to study the 
time dependence Wi{t). Strictly speaking, this formula 
is valid near the center of the energy spectrum, otherwise 
one should take into account distortion effects near the 
edges of the energy spectrum. 

Due to the normalization conditions, / P{E)dE = 1 
and / E'^P[E)dE = A|;, we have the following relations 
El, 
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and 
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where $(z) is the error function. 

The return probability W{t) corresponding to the 
strength function (p7|) is defined by the integral. 



A{t) = B dE 



exp 
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(20) 




FIG. 1. Schematic time dependence W{t) for f = 0.5, 
Ae = 1.2; the dependence W{t) = exp{ — A%t^) changes into 
W{t) = exp{-rt) at the point tc = Tp/A% ^ 0.17. Simple 
exponential dependence is shown by the dotted line. 



This qualitative result is supported by a more detailed 
consideration. In many-body systems the dependence 
F(i?) can be approximated by the Gaussian function 
since the density of final states pf{Ef) in Eq.(|^) is typi- 
cally close to the Gaussian Q]. This provides us with an 
approximate formula for the strength function (tested by 
the numerical calculations 



(E^Eo) 



The time dependencies of A{t) and W{t) for small time 
are given by Eqs.(^3[|lJ). If F ^ cr, the region of the ap- 
plicability of these equations is very narrow. Indeed, in 
this case A|; « and the condition t <^tc<^ ~ ^ 

results in the relation A\t^ <C 1 . The absolute value of 
the amplitude A{t) in this case is given by the series in 
the parameter {atY = {i/icY i 
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For large time, t ^ tc = 1/cr, the calculation of the 
integral in Eq.(20) gives 



Wit) « exp 



ill 

ttAI 



Ft 



for the return probability. Here the correction -^-^ ~ 
— 2£= is small. Indeed, the strength function in this case 
is close to the Lorentzian which gives a simple depen- 
dence W{t) — exp(— Ft) for the decay. 

Another limit case of a strong interaction, F ^ A^ 
(or, the same, F ^ tr), is more delicate. In this case the 
strength function is close to the Gaussian with Ae ~ cr 
and tc is large, tc ^ ^ ^. The leading dependence of 
W{t) in this case is the Gaussian, W{t) ~ exp{—A'^t'^) 
. Only for a long time t ^ ^ it becomes the simple 
exponential function. 



It is important to stress that even for a large time 
the return probability W{t) has large correction factor 

exp (AE)2 ) ' addition to the standard decay law 
exp(-ri). 

Due to a finite number of particles, there are spe- 
cific important features in the dynamics of wave packets, 
namely, the damped oscillations and the break of the de- 
cay for Wi{t) The number of basis components |/) 
within the energy shell \Eq — Ef \ < min(r, cr) = A is 
finite. Therefore, the decay stops if Wj is clo se to the 
equilibrium value defined as Woo = Wi{t ^ oo) 3N~^^. 
Here Npc is the number of principal components in an 
eigenstate, Npc ^ A/D , where D = is the mean en- 
ergy interval between all many-body levels. Note that 
the value of Woo is still at least 3 times larger than 
Wf = N~c^ for any other component / , W,;(cxd) > 

Wf {oo) (see, details in |l5|]). 

The equilibrium occurs because the average decay fiux 
is equal to the average return flux. However, the re- 
turn flux leads to the damped oscillations of Wi{t) and 
to the oscillations of a current number of the princi- 
pal components Npc{t). These oscillations arise because 
the decay flux "reflects" from the edges of the energy 
shell, when all components within this shell are popu- 
lated. Period of these oscillations is about ric/A where 
A is the inverse decay time, and ric is the number of 
"classes" in the Hilbert space. This number can be de- 
fined as the number of interaction steps in the perturba- 
tive chain H^^a^Ha^a^ . . . Ha^a„^, needed to populate all 
basis states within the energy shell. For example, in the 
TBRI model with 6 particles and 12 orbitals, the num- 
ber of steps is ric = 3 since each two-body interaction 
Hik moves only two particles to new orbitals. 

In conclusion, we have studied generic features of the 
return probability W(t) for a system to be found in an ini- 
tially excited many-body state. Due to the two-body in- 
teraction between Fermi-particles, the wave packet in the 
energy representation spreads over all basis states within 
the energy shell. The dependence W{t) for small time 
is determined by a ballistic spread of the packet and is 
given by the expression For large time, the decrease 
of W{t) is determined by the form of the strength func- 
tion P{E) . We have analyzed the behavior of W{t) by 
making use of the analytical expression for P{E), which 
is obtained for any strength of random two-body interac- 
tion between finite number of interacting Fermi-particles. 

We have shown that for the Breit-Wigner form of P{E) 
(relatively weak interaction) the decay of W{t) on a large 
time scale has the conventional exponential dependence, 
W{t) ~ exp(— Ft). On the other hand, for the Gaus- 
sian form of P{E) (strong interaction) the time depen- 
dence W{t) turns out to be of very specific. Namely, the 
leading term gives the quadratic exponential dependence, 
W{t) ^ exp{—A^t^), and only for a very large time the 



conventional exponential linear dependence formally re- 
covers. However, in this case an additional prefactor C 
appears before the exponent, which turns out to be ex- 
ponentially large, thus resulting in a strong modification 
of the standard exponential estimate for W{t). 
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